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Abstract 

We explore the sensitivity of W and Z boson production in hadronic collisions to uncertainties 
in parton distribution functions (PDFs) at large x arising from uncertainties in nuclear corrections 
when using deuterium data in global QCD fits. The W and Z differential cross sections show in- 



creasing influence of nuclear corrections at high boson rapidities, particularly for the d quark, which 
7^ , ' is diluted somewhat in the decay lepton rapidity distributions. The effects of PDF uncertainties 

on heavy W' and Z' bosons beyond the Standard Model become progressively more important 
for larger boson masses or rapidities, both in pp collisions at the LHC and in pp scattering at the 
Tevatron. 
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I. INTRODUCTION 



The discovery and determination of properties of new particles beyond the Standard 
Model at high-energy colliders depends on accurate knowledge of parton distribution func- 
tions (PDFs) of the hadrons involved in the coUisions. With the Large Hadron Collider 
(LHC) at CERN now taking data at unprecedentedly high energies, the effort to control back- 
grounds in searches for the Higgs boson and other putative particles is taking on paramount 
urgency. Recently the dependence of the Higgs boson cross section for the dominant gluon- 
gluon fusion channel on PDFs has been the cause of some debate [1-3], highlighting the 
need for a careful determination of strong interaction inputs such as gluon distributions, the 
strong coupling constant cKs, and higher-order radiative corrections. 

At lower energies, the importance of PDF uncertainties has also been discussed recently 
in fixed-target experiments, particularly for the d quark distribution in the region of large 
parton momentum fractions x [x 0.5) [4, 5]. Because both proton and deuterium deep- 
inelastic scattering (DIS) data are required to constrain the d quark PDF, uncertainties 
in the nuclear corrections in the deuteron at large x translate into significant and growing 
uncertainties on the d/u ratio as x — > 1. Through evolution, this can impact cross section 
calculations at smaller x and larger Q"^ [6] , especially in regions where the rapidity is large. 

In addition to Higgs boson cross sections, other processes studied at the LHC or the 
Tevatron at Fermilab that may be sensitive to PDF uncertainties include the production of 
heavy W and Z' bosons associated with additional SU(2)xU(l) gauge groups. These are 
predicted in various extensions of the Standard Model, such as the SO(IO) and Eg grand 
unified theories, or supersymmetric models, some with W and Z' boson masses at the TeV 
scale (for reviews see, e.g., Refs. [7-10]). Their production cross sections at large rapidities 
will involve products of PDFs evaluated with one value of x small and the other large, 
thereby exposing them to uncertainties in PDFs at large x, particularly near the kinematic 
limits. 

In this paper we explore the sensitivity of the weak boson production cross sections to 
uncertainties in PDFs at large x. For our numerical estimates we use the PDFs from the 
recent CTEQ- Jefferson Lab (CJ) next-to-leading order (NLO) global analysis [5] of proton 

and deuteron data, which quantified the model dependence of the nuclear corrections in the 
deuteron and the resulting effects on the PDFs. In Sec. II we briefly review the CJ analysis 
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and the origin of the PDF uncertainties, before examining their impact in Sec. Ill on the 
physical and Z boson cross sections, and W and lepton charge asymmetries. The effects 
of large-x PDF uncertainties on production rates of heavy W and Z' bosons in pp and pp 
coUisions at the LHC and Tcvatron, respectively, are studied in Sec. IV as a function of the 
boson mass, and limits placed on the accuracy with which cross sections for bosons of a given 
mass can currently be determined. (Even though data taking has now been completed at the 
Tevatron, considerable quantities of data remain to be analyzed.) Finally, some concluding 
remarks are made in Sec. V. 

II. PDF UNCERTAINTIES AT LARGE X 

The CJ analysis [5] was a dedicated global NLO fit of proton and deuteron DIS and other 
high-energy scattering data, which critically examined the effects on PDFs of nuclear cor- 
rections in the deuteron F2 structure function. Nuclear corrections were estimated using a 
smearing function computed within the weak binding approximation [11, 12], taking into ac- 
count nuclear binding and Fermi motion, as well as a range of models describing the possible 
modification of the nucleon structure function off-shell [11, 13, 14]. Several deuteron wave 
functions were considered, based on high-precision nucleon-nucleon potentials, including the 
nonrelativistic CD-Bonn [15] and AV18 [16] wave functions, and the relativistic WJC-1 and 
WJC-2 wave functions [17], as well as the older Paris [18] wave function for reference. The 
off-shell corrections were estimated from a relativistic quark spectator model [14] , and from 
a phenomenological model proposed by Kulagin & Petti [11] but modified for the specific 
case of the deuteron (see Ref . [5] for details) . 

Combinations of deuteron wave functions and off-shell models giving the smallest and 
largest nuclear effects were identified, and used to define the range of the nuclear corrections 
from the minimum (WJC-1 wave function and no off-shell corrections) to the maximum 
(CD-Bonn wave function and largest off-shell corrections) nuclear corrections. The central 
values, which are used as a reference, were obtained using the AVIS wave function and an 
intermediate off-shell correction. The resulting fitted d/u quark distribution ratio is shown in 
Fig. 1 for the full range of nuclear uncertainties determined in Ref. [5] (sec also Rcfs. [19, 20]). 
While the u quark distribution is relatively well constrained by proton structure function 
data for all values of x, the d quark PDF has large uncertainties beyond x ~ 0.5. Significant 
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FIG. 1: Ratio of d to n quark distributions for the CJ PDFs [5] at = 10 GeV^. The shaded band 
ihustrates the uncertainty range between maximum (blue dashed) and minimum (red dot-dashed) 
nuclear corrections in the deuteron. 

uncertainties appear also in the gluon PDF in the region of x not directly constrained by 
data, as well as in the u and d distributions due to correlations with the d quark, induced 
in the global fits by moderate-x jet and dilepton production data, respectively. 

Several experiments to directly measure (i/w up to x ~ 0.8 are planned at the 12 GeV 
energy upgraded Jefferson Lab in the near future [21-23], which it is hoped will reduce the 
uncertainties significantly. These include the "MARATHON" [22] experiment, which aims 
to extract jF^ from a measurement of the F^ structure functions of tritium and ^He with 
cancellation of nuclear effects to the ^ 1% level [24], as well as the "BoNuS" experiment 
[21], which minimizes nuclear corrections in semi-inclusive DIS from deuterium by tagging 
slow, backward protons that effectively guarantee scattering from a nearly-free neutron. In 
addition, parity-violating DIS on a hydrogen target [23] will yield a new combination of u 
and d PDFs at large x, free of any nuclear corrections (see Ref. [5] for further details). In the 
meantime, however, it is important to establish the limitations that the current uncertainties 
on the d quark PDF at large x place on the calculation of observables which may be sensitive 
to these. In the following we shall illustrate the impact of PDF uncertainties at large x arising 
specifically from the dependence on the model of nuclear corrections in the deuteron. The 
choice of CJ PDFs for this purpose is merely for convenience, as these are the only PDFs 
available that explicitly quantify the nuclear model dependence. 
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III. W AND Z BOSON PRODUCTION 



In this section we discuss the effects of PDF uncertainties at large x on the W and 
Z boson cross sections, and possible constraints on these obtained from measurements at 
large rapidities. Earlier studies probing the sensitivity of weak boson production to PDF 
uncertainties were explored in Refs. [25-33]. The discussion here is not meant to provide 
an exhaustive account of detailed aspects of W boson production, but simply highlight the 
fact that nuclear corrections in deuterium are an important source of PDF uncertainty at 
large x that has not been addressed in earlier analyses. To begin with we shall review the 
general formulas for the inclusive weak boson production cross sections in hadronic collisions 
relevant to current collider experiments. 

A. Cross sections 

Hadron-hadron collisions involve at least two interacting partons, one from the hadron 
"beam" and one from the "target", with momentum fractions xi and X2, respectively. At 
fixed center of mass energy y/s and boson rapidity 



where E and Pz are the boson energy and longitudinal momentum in the hadron center 
of mass frame, the par ton momentum fractions are given (at leading order in the strong 
coupling constant) by 



where M is the mass of the produced boson. The absolute value of the rapidity thus ranges 
from up to \y\raax — log(v^/M). For inclusive production in pp or pp collisions, for 
example, the cross sections (to leading order and neglecting heavy quarks) are given by [34] 




(1) 



Xl,2 = 




(2) 



da 
dy 



{pp W+X) 



2nGF 
3\/2 



X1X2 (cos^^c [u{xi)d{x2) + d{xi)u{x2)] 
+ sin^^c [u{xi)s{x2) + s{xi)u{x2)]) , 

X1X2 (cos^ 9c [u{xi)d{x2) + d{xi)u{x2)] 
+ sin^ 6c [u(xi)s(x2) + s(xi)u(x2)]) , 



(3a) 
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where is the Fermi constant, and 9c is the Cabibbo mixing angle. The W~ differential 
cross sections are similar to those in Eqs. (3), but with quark PDFs replaced by the corre- 
sponding antiquark PDFs. Consequently the cross sections in pp collisions are related 
by {daw+/dy){y) = {day/- / dy){—y), and hence are equivalent when integrated over rapidity. 
For pp collisions the individual and W~ cross sections are symmetric with respect to 
y ^ —y, but otherwise unrelated. 

Similarly, the leading order, light quark cross sections for Z boson production in pp or 
pp collisions are given by [34] 

(pp ^ ZX) = ^ Yl [(9vf + i9lf] ^1 ^2 + q{^i)q{x2)) , (4a) 



dy 3V2 
da , 2txG 



1 



{pp ^ZX)^ [(9vf + ig^f] ^1 ^2 {q{x,)q{x2) + q{^i)q{x2)), (4b) 

where Qy — — 2eq sin^ 9w and g\ — are the vector and axial- vector couplings of the 
Z boson to quark q [35], with and tl the electromagnetic charge and weak isospin of 
the quark, respectively, and 6w the weak mixing angle. The symmetry properties of the 
differential Z cross sections are such that {daz / dy){y) = {daz / dy){—y) for both pp and pp 
collisions. Note that in the conventions of Ref. [34] the couplings gy j^ are two times smaller 
than the standard ones in Ref. [35]. In the convention used here the couplings {gyY -\- {q'a}'^ 
in Eqs. (4) are equal to 5/18 + A(l + A)/9 and 13/36 + A(l + A/4)/9 for u and d quarks, 
respectively, where A = 1 — 4sin^^iy. Since A 0, the effective strengths of the Z boson 
couplings to u and d quarks are therefore similar. 

While the expressions in Eqs. (3) and (4) are given at leading order, in practice we 
compute all cross sections at NLO, including heavy quarks. The leading order expressions 
give the dominant contributions, however, and are instructive in clearly illustrating that in 
pp collisions, for instance, the cross section at large rapidity is mostly dependent on the 
u quark, while the W~ cross section depends mostly on the d. Since the nuclear corrections 
discussed in Sec. II induce the greatest uncertainty into the large- a; d quark PDF, one can 
immediately deduce that production in pp scattering will be most affected by these 
uncertainties, while production will be relatively inert. For pp collisions, the large-x 
PDF uncertainties will affect W~ cross sections at large positive rapidities, or equivalently 
W'^ cross sections at negative rapidities. 

In the following sections we will compute the W and Z boson cross sections numerically 
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to study their sensitivity to PDF uncertainties at large x. All calculations will be for 
pp collisions at the LHC with ^/s = 7 TeV, and for pp collisions at the Tevatron with 
y/s = 1.96 TeV. 



B. Z bosons 

The sensitivity of the differential Z boson cross section to the different PDF behaviors 
at large x is illustrated in Fig. 2 as a function of the Z boson rapidity yz, for LHC and 
Tevatron kinematics. The cross sections are computed from the CJ PDFs [5] with minimal 
and maximal nuclear corrections, relative to a reference cross section computed from the 
central PDFs as described in Sec. II. 
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FIG. 2: Differential Z boson cross section as a function of the Z rapidity yz, computed from CJ 
PDFs with maximum (blue dashed) and minimum (red dot-dashed) nuclear corrections, relative to 
the reference cross section (T^(ref) calculated using the central CJ PDF set [5]. The cross sections 
are computed for pp collisions at the LHC with ^/s = 7 TeV (left) and for pp collisions at the 
Tevatron with ^/s = 1.96 TeV (right). 



The behavior of the cross section ratios is qualitatively similar at both the LHC and 
the Tevatron, with the main difference being the range of rapidities accessible from the 
respective available energies ^/s. At low rapidities the cross sections are relatively insensitive 
to uncertainties in the large-x PDFs, with differences of ^ 1% for yz ^ 3 at the LHC and 
yz 2 at the Tevatron. At larger rapidities, however, there is far greater sensitivity to the 
large-x behavior, particularly of the d quark, leading to ~ 15% uncertainty in the differential 
cross section for j/^ = 4 at the LHC, and for yz = 2.8 at the Tevatron, which correspond to 



parton fractions of x 0.7. As one approaches the kinematical thresholds of yz,msx ~ 4.3 
at the LHC and yz,max ~ 3.1 at the Tevatron, these uncertainties increase dramatically, as 
would be expected from the x — > 1 behavior of the d/u ratio in Fig. 1. 

C. W cross sections and asymmetries 

The behavior of the W boson differential cross sections as a function of the W rapidity 
Uw is qualitatively similar to those for the Z, but with some important differences, as Fig. 3 
illustrates for W'^ and W~ production at LHC and Tevatron kinematics. Again, there 
is very little dependence on the large-x PDF uncertainties at low rapidity, but increasing 
sensitivity as the rapidity approaches its kinematic upper limit of yw,maK ~ 4.5 at the LHC 
and yz,max ~ 3.2 at the Tevatron. 

For bosons, the cross section maintains relatively little dependence on the large-x 
nuclear corrections over the entire rapidity range, barely reaching 4% difference at yw ~ 4 
at the LHC, or yw 3 at the Tevatron. In contrast, the W~ cross section shows an even 
stronger dependence on nuclear corrections than the Z cross section in Fig. 2, deviating 
significantly from unity for yw ^ 3 for the LHC and yw <^ 2 for the Tevatron, and reaching 
upwards of 40% deviation at yw ~ 4 and 3 for LHC and Tevatron kinematics, respectively. 
The greater sensitivity of the W~ production cross section compared with the can 
be understood from the dominance of the latter by the u quark PDF at large x, which 
is relatively insensitive to the nuclear correction uncertainties. The enhancement of the 
W~ cross section at large yw for the CJ PDFs with maximum nuclear corrections, and 
corresponding suppression of the CJ PDFs with minimum nuclear corrections, relative to 
the central CJ fits essentially follows the trend of the d quark PDF in Fig. 1. The slight 
enhancement of the cross section at large yw for the CJ PDFs with minimum nuclear 
corrections reflects the anticorrelation of the u quark PDF with respect to the d observed 
in Ref. [5]. 

Taking differences and sums of the and W~ cross sections, one can construct the W 
boson asymmetry, 

. _ <^w+{y) -(^w-{y) 

— ;r^r- 1-^, [o) 

crw+{y) + crw-{y) 

where aw±iy) = daw±/dy. The asymmetry is shown in Fig. 4 versus the W rapidity at 
the LHC and Tevatron for the CJ PDFs with maximum and minimum nuclear corrections. 
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FIG. 3: Differential (top) and (bottom) boson cross sections as a function of the W 
rapidity yw, computed from CJ PDFs with maximum (blue dashed) and minimum (red dot-dashed) 
nuclear corrections, relative to the reference cross sections cr^^/i (ref) calculated using the central 
CJ PDF set [5]. The cross sections are computed for pp collisions at the LHC with ^/s = 7 TeV 
(left) and for pp collisions at the Tevatron with ^/s = 1.96 TeV (right). 

A clear deviation between the two PDF fits becomes visible at yw ^3.5 for the LHC and 
yw ~ 2 for the Tevatron, corresponding to one of the partons carrying momentum fractions 
X ~ 0.4 and x ~ 0.35, respectively. Data on W boson asymmetries may therefore provide 
constraints on the d/u quark distribution ratio already at these moderate values of x. In 
particular, comparison with the CDF data [36] in Fig. 4 illustrates a preference for larger 
Aw values at yw ~ 2, which corresponds to xi ^ X2- As observed in Ref. [5], since the 
asymmetry at large yw can be approximated by 

d{x2)/u{x2) - d{xi)/u{xi) 



A 



w 



[Xi > X2] 



(6) 



d{x2)/u{x2) + d{xi)/u{xi) ' 
this would suggest a smaller d/u ratio at large Xi, as would arise for the CJ PDFs with 
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FIG. 4: W boson asymmetry Ay^ as a function of the W rapidity yvK a-t LHC (left) and Tevatron 
(right) kinematics, computed from CJ PDFs [5] with minimum (upper blue sohd) and maximum 
(lower blue solid) nuclear corrections. For comparison, the asymmetries using the ABKM [29] (red 
dot-dashed), CTIO [30] (green dashed) and MSTW [31] (black dotted) PDF sets are also shown. 



minimum nuclear corrections. 

For comparison, the asymmetries calculated from the ABKM [29], CTIO [30] and MSTW 
[31] PDF sets are also shown in Fig. 4. Each parametrization shows good agreement with 
the Tevatron CDF data, with the exception of the ABKM fit, which overestimates the 
asymmetries at intermediate rapidities, yw ~ 1 — 2. This may be due to the W asymmetry 
data not being fitted directly in the ABKM analysis. At large rapidity the spread in the 
various PDF sets is comparable to the difference between the CJ PDFs with minimal and 
maximal nuclear corrections. However, we stress that the origin of the differences between 
the PDF sets is unrelated to the difference between the two CJ PDF sets. Had the various 
non-CJ PDFs sets included nuclear uncertainties in their analysis, each one would have 
a corresponding nuclear uncertainty band similar to the one in Fig. 4, and the combined 
spread between them would subsequently be significantly larger. 



D. Lepton Asymmetries 

Experimentally, measurement of W bosons asymmetries requires reconstruction of the W 
boson distributions from their leptonic decays, — > l^vi and W~ — t- l~vi^ with / = e or 
/i. On the other hand, lepton charge asymmetries can be constructed directly from the 
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decay products and studied as a function of the lepton pseudorapidity t], defined as 



1 



V 



In 



Ifcl + k. 



— In tan ■ 



e 



(7) 



2 V I I A^2: J 2 

where k is the charged lepton momentum, and 6 is the angle between the lepton momentum 
and the beam axis in the center of mass frame. The lepton asymmetry is then given by 



(8) 



where ai± {rj) = dai± / dr] is the differential cross section for the production and leptonic decay 
of the W^. 

Lepton asymmetry data from the DO Collaboration [37] at Fermilab are shown in Fig. 5 
as a function of the pseudorapidity up to ?7 ~ 2.5. The data are integrated over lepton 
transverse momenta pt > 25 GeV, and compared with asymmetries computed from the CJ 
PDFs with maximum and minimum nuclear corrections using the MCFM (Monte Carlo for 
FeMtobarn processes) program [38]. Good agreement is obtained between the calculated 
asymmetry and data, although little sensitivity is evident to the large-x nuclear uncertainty 
in the PDFs observed in Fig. 4 until t] ^ 3. 
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FIG. 5: Lepton charge asymmetry A^j as a function of lepton pseudorapidity rj at the LHC (left) 
and Tevatron (right), computed from CJ PDFs [5] with maximum (blue dotted) and minimum 
(red solid) nuclear corrections. The calculations are compared with CMS W — )• ei^e (red squares) 
and W — )• fiv^ (blue triangles) data for a lepton transverse momentum cut px > 25 GeV [40], and 
preliminary LHCb data (black circles) for pT > 20 GeV [41], as well as with DO data from the 
Tevatron pr > 25 GeV [37]. 
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Similar behavior is found for the new lepton charge asymmetry data from the LHC. Here 
lepton asymmetries have been measured by the ATLAS [39] and CMS [40] Collaborations for 
pseudorapidities |?7| ^ 2, and preliminary results from the LHCb Collaboration [41] extend 
the coverage to 2 ^ 77 < 4.5 [42]. The agreement of the CJ PDFs with the LHC data is 
good over the entire range of rj, as Fig. 5 illustrates for CMS W ^ ev and W ^ fiu data 
integrated over lepton transverse momenta pr > 25 GeV, and LHCb data with pr > 20 GeV. 
The dependence on the large-x behavior of PDFs, however, becomes visible only for 77 ^ 4. 

The limited sensitivity of the lepton asymmetries to large-x PDFs is not surprising, given 
that the lepton asymmetry is computed by convoluting the W boson cross sections with the 
W boson decay distributions, which dilutes the sensitivity to regions where the PDFs are 
small. Although the lepton asymmetry data are clearly valuable for constraining global PDF 
fits in general, greater sensitivity to the large-x behavior of the d/u ratio may be possible 
through the reconstruction of the W boson asymmetries themselves. Note, however, that 
the reconstruction of W boson asymmetries is limited by theoretical uncertainties such as 
the modeling of the pr distributions and higher-order resummation corrections (which also 
affect the lepton asymmetries), as well as the choice of PDFs used to compute the event 
reweighting coefficients in the reconstruction. 

IV. HEAVY W AND Z' BOSONS 

The production rate of any new heavy boson beyond the Standard Model will naturally 
depend on its internal properties such as the spin. Many possibilities have been canvassed 
for how such heavy bosons can arise [9] , including as scalar excitations in i?-parity violating 
supersymmetry [43], spin-1 Kaluza-Klein excitations of Standard Model gauge bosons in 
the presence of extra dimensions [44], or as spin-2 excitations of the graviton [45]. On the 
other hand, if the new bosons are associated with extensions of the Standard Model gauge 
group, their interactions with fermions will resemble those of the W and Z bosons of the 
electroweak theory, with different masses and couplings. 

In this section we explore this latter possibility, and in particular the sensitivity of the 
production cross sections to uncertainties in PDFs at large x. From Eq. (2) one can see 
that increasing the mass will directly increase the relevant x values, so that higher mass 
bosons will more readily sample the high-x region where the nuclear uncertainties are more 
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prominent. In the calculations discussed here we shall assume that the putative W and Z' 
bosons have the same properties as the Standard Model W and Z bosons, except for their 
larger masses. The cross sections will of course decrease rapidly with increasing boson mass, 
so that the effects of the large- a; PDF uncertainties will become more significant as the mass 
increases. Of course the details of the predictions will change in more sophisticated models 
in which the W and Z' couplings are different from those in the Standard Model; however, 
the simplified scenario considered here is sufficient to illustrate the possible impact of large-x 
PDF uncertainties on new physics searches. 

Currently the exclusion limit on the W mass from the Tevatron, for couplings similar to 
Standard Model couplings, is Myyi > 1.12 TeV at the 95% confidence level with an integrated 
luminosity of 5.3 fb~^ [46]. For the neutral Z' boson the limits vary between Mz' ^ 800 GeV 
and Mz' ^ 1 TeV [47], depending on the Standard Model extension considered [35]. The 
latest results from the LHC place the limits for the W mass at Mw' > 2.15 TeV [48] and 
for the Z' mass at Mz' > 1.83 TeV [49] in the Sequential Standard Model, with the same 
couplings to fermions as for the W and Z. A similar constraint on the Z' mass in grand 
unified theories is also obtained from measurements of atomic parity violation in ^^^Cs [50]. 

A. Rapidity distributions 

The differential cross sections for heavy W and Z' bosons will generally have similar 
behavior as a function of rapidity to those of the physical boson cross sections in Figs. 2 
and 3, except for a smaller rapidity range, with the large-x region emphasized more strongly 
for increasing boson mass. According to the study by Erler et al. [51], pp collisions at 
the LHC with ^/s = 1 TeV will be sensitive to Z' masses up to Mz> ^2.1-2.7 TcV for 
luminosities between 30 fb^^ and 300 fb~^, corresponding to the so-called "low-luminosity" 
and "high-luminosity" LHC scenarios, respectively. (With ^/s — 14 TeV the mass limits 
would vary between Mz' ~ 3.6 and 4.6 TeV.) For pp collisions at the Tevatron, an energy 
of y/s — 2 TeV with a luminosity of 10 fb~^ would be expected to allow sensitivity to Z' 
masses up to f« 1 TeV [51]. In this section we therefore consider Z' (and W) masses up to 
1 TeV and 3 TeV for Tevatron and LHC kinematics, respectively. 

From Eq. (2), larger boson masses naturally restrict the kinematically accessible range of 
rapidities, so that at the LHC, for example, a 1 TeV (3 TeV) Z' boson can be produced at 
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FIG. 6: Differential Z' (top), W'^ (center) and W'~ (bottom) cross sections as a function of 
the rapidity, computed from CJ PDFs with maximum and minimum nuclear corrections, relative 
to the reference cross sections (T^/^iy'(i"ef) calculated using the central CJ PDF set [5]. The LHC 
cross sections (left) are computed for boson masses Mz'^v = 1 TeV (red dot-dashed), 2 TeV 
(short-dashed) and 3 TeV (long-dashed) with = 7 TeV, while the Tevatron ratios (right) are 
shown for Mz\w' = 0-5 TeV (red dot-dashed), 0.75 TeV (short-dashed) and 1 TeV (long-dashed) 
with ^ = 1.96 TeV. 
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a maximum rapidity of |yz'|max = 2.0 (0.9), compared with |yz|max = 4.3 for the Standard 
Model Z boson. The kinematic reach and sensitivity to large-x PDFs is illustrated in 
Fig. 6 (top) for the differential Z' cross section ratio as a function of yz'-, for masses Mz' = 1, 
2 and 3 TeV at the LHC, and M^/ = 0.5, 0.75 and 1 TeV at the Tevatron. As the rapidities 
approach their kinematical thresholds for a given M^/, the uncertainty in the differential 
cross sections increases significantly, reaching about 30 — 40% of the central CJ value for pp 
coUisions at the LHC, and 15 — 20% for pp coUisions at the Tevatron. 

Because the couplings of the Z' to u and d quarks are assumed to be similar (see 
Sec. Ill A), the Z' cross section at the LHC depends on both the combinations u{xi)u{x2) + 
u{xi)u{x2) and d{xi)d{x2) + d{xi)d{x2)- However, since the d/u and u/u ratios arc <^ 1 at 
large x values, which are preferentially sampled for large M^/, the contributions of d quarks 
are suppressed relative to u quarks. Consequently the Z' ratios in Fig. 6 are only mildly 
affected by uncertainties in quark PDFs at large x. The Z' production cross section in pp 
collisions at the Tevatron, on the other hand, is determined predominantly by the product 
u{xi)u{x2)i which is well constrained and independent of the nuclear model for all xi,2, and 
therefore has an even smaller uncertainty. 

For the W differential cross sections, the behavior as a function of yw' is qualitatively 
different for W'^ and W'~ production, shown in Fig. 6 (center) and (bottom), respectively, 
with the latter displaying dramatically greater sensitivity to large-x PDF uncertainties. 
This is clear from Eqs. (3), where for pp coUisions the dominant contribution to the W'^ 
cross section depends on the products u{xi)d{x2) and d{xi)u{x2). While the u quark PDF 
is insensitive to the nuclear corrections, the d distribution varies considerably with the 
nuclear model, especially at larger values of x. At high rapidity the d/u ratio is small, 
and the cross section is determined by the u PDF with xi large and the d PDF with 
X2 small, both of which are well constrained. For y-^i — > 0, on the other hand, one has 
X2^ 0.14 for Mw = 1 TeV and xi ^ X2 0.42 for Mw = 3 TeV, at which the d 
PDF has significantly greater uncertainty than the u, yielding up to ^ 25% uncertainties 
in the cross section. Similarly for pp collisions at the Tevatron, the W'^ cross section is 
dominated by the combination u{xi)d{x2)., which for Xi ^ X2 at high rapidity is relatively 
well constrained. At central rapidity, with Xi = X2^ 0.25 (0.5) for M^y = 0.5 (1) TeV, the 
uncertainties in the cross section remain within the ~ 10% level. 

In contrast, the W'~ cross section at high rapidity is determined by the products 
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d{xi)u{x2) and u{xi)d{x2) in pp collisions, and d{xi)u{x2) and u{xi)d{x2) in pp collisions. 
Consequently, uncertainties in the cross sections at large rapidity, both at the LHC and the 
Tevatron, arise mainly from the d quark at large x, and exceed 100% as the kinematic limit 
in yyi/ is approached. Qualitatively, the growing uncertainty of the W'^ cross section with 
increasing rapidity resembles the W~ cross section ratio at large y^r in Fig. 3. At central 
rapidity, the uncertainty in the pp cross section at the LHC is of the order 10%, arising 
mainly from the uncertainty in the u distribution. At the Tevatron, the pp cross section 
is well constrained for small boson mass, but for Mw' > 0.75 TeV, with xi — X2 ^ 0.3, 
becomes increasingly sensitive to uncertainties in the d quark PDF, reaching about 20% for 
Mw' = 1 TeV. 

B. Integrated cross sections 

Integrating over all rapidities, the resulting total Z' cross section computed from CJ 
PDFs with minimum and maximum nuclear corrections is shown in Fig. 7 as a function of 
the Z' mass. Relatively little dependence on the PDFs is observed, with effects of ^ 3% 
observed for Mz' < 3 TeV at the LHC and Mz' < 1 TeV at the Tevatron. This is not 
surprising given that total cross sections are dominated by contributions from low values of 
yz', where the PDF uncertainties are generally smaller than at high values of yz', at which 
the contributions are suppressed by the steeply falling PDFs as x ^ 1. At larger Mz' values 
the uncertainties generally increase, but are subject to greater fluctuations in the antiquark 
distributions at high x, and hence are less reliable. 

The integrated cross sections for W bosons in Fig. 8 show somewhat greater sensitivity 
to large-a; PDFs as a function of M^^r. For W'~^ bosons produced in pp collisions at the LHC 
the uncertainties increase from ^ 5% for Mz' = 1 TcV to ~ 20% for Mz' = 3 TeV. This 
behavior stems directly from the increasing uncertainty in the d antiquark PDF at large x 
apparent in the W'^ rapidity distribution at low yw' in Fig. 6. For W'~ boson production 
in pp scattering, the uncertainties in the total cross section are smaller than for the W'^ at 
low Mw', remaining < 2% for Mw' < 2 TeV, but increase to fa 10% at Mw' = 3 TeV due 
to the uncertainty in the u quark. The stronger dependence on the behavior of the d quark 
PDF at large x apparent in the W'^ differential cross section at high rapidity in Fig. 6 is 
mostly washed out in the integrated cross section. 
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Because the W cross sections in pp collisions at the Tevatron are determined by the 
products u{xi)d{x2) and d{xi)u{x2) for W'^ and W'~ , respectively, integrating over rapidity 
samples all accessible values of Xi and X2, so that the total W'~^ and W'~ cross sections are 
equivalent. The dependence of the integrated W cross sections on PDFs essentially follows 
the d quark distribution. For masses M^y/ ^ 0.5 TeV there is little sensitivity to the large- 
X behavior of the PDFs, with ^ 5% uncertainty in the cross section ratio, but increasing 
dependence at larger Mw', with ^ 30% uncertainty at Mw' = 1 TeV. 

While the sensitivity of the W and Z' cross sections to the large-x behavior of PDFs 
increases with increasing W and Z' masses, the absolute values of the cross sections naturally 
fall with increasing masses, some 3 orders of magnitude from 100 GeV to 3 TeV. This is 
illustrated in Fig. 9, where the ratio of the integrated W'~^ + W'~ cross sections computed 
from CJ PDFs with minimum and maximum nuclear corrections, relative to the cross section 
with the central CJ PDFs, is plotted versus the integrated Z' cross section. Here the ratio 
of the W to Z' masses is kept constant in order to study the effect of the increasing W, Z' 
mass. For larger boson masses the impact of the large-x PDF uncertainties clearly increases, 
reflecting the trend observed in Figs. 7 and 8. Note that because the integrated W'^ cross 
section is generally larger than the W'~ cross section (because of the larger u distribution 
compared with the d), the cxvy /avy (ref) ratio in Fig. 9 generally follows the ratio of the W'~^ 
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FIG. 7: Integrated Z' boson cross section from Fig. 6 as a function of the Z' mass, computed 
from CJ PDFs with minimum (red dot-dashed) and maximum (blue dashed) nuclear corrections, 
relative to the reference cross section az' (ref ) calculated using the central C J PDF set [5] , for LHC 
(left) and Tevatron (right) kinematics. 
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cross sections in Fig. 8 for increasing boson mass. 



V. CONCLUSION 



In this paper we have explored the sensitivity of weak boson production in hadronic 
collisions to parton distributions at large values of x. At present there are large uncertainties 
in the d quark distribution, particularly above x ~ 0.5, arising from the model dependence 
of nuclear corrections used when analyzing deuteron DIS data in global PDF fits, which 
can impact cross section measurements at large rapidities. The PDF uncertainties can also 
affect production cross sections of heavy W and Z' bosons beyond the Standard Model at 
central rapidities. 

Using PDFs extracted from the recent CJ global fit [5], we find increasing sensitivity to 
the large- a; region for Z boson production in pp collisions at the Tevatron for Z rapidities 
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FIG. 9: Ratio of W' boson cross sections aw' /crw'{^^^)i where aw' = o'wi+ + cr\y/-, versus the Z' 
boson cross section az' for varying Z' boson masses, indicated by the numbers (in TeV) next to the 
points (the points corresponding to the physical W and Z cross sections are labeled by ^W,Z"). 
The cross sections are computed from CJ PDFs with minimum (filled red circle) and maximum 
(open blue circle) nuclear corrections, relative to the reference cross section calculated from the 
central CJ PDFs [5]. 

yz ^ 2, and in pp collisions at the LHC for ~ 3. Precision measurements of Z boson cross 
sections at these rapidities, at the Tevatron and particularly at the LHC with the LHCb 
experiment, will be required to impact global fits of the d quark and constrain nuclear model 
uncertainties. 

For charged weak bosons, the cross sections are mostly independent of PDF uncer- 
tainties due to their preferential coupling to u quarks, whereas the W~ cross sections display 
strong dependence on the d quark uncertainties for W rapidities yw ~ 1-5 at Tevatron and 
yw ~ 3 at LHC kinematics. Measurements of W boson charge asymmetries at large ra- 
pidities, such as those from the CDF Collaboration at Fermilab [36], thus provide strong 
constraints on the behavior of the d/u ratio at large x, although such measurements are 
very challenging given the low rates expected in the relevant regions of kinematics. Direct 
reconstruction of W boson asymmetries in the LHCb experiment for ^/s = 7 TeV would 
also be extremely valuable in providing access to PDFs at x ~ 1. We have also compared 
charged lepton asymmetries with data from DO [37], and from the CMS [40] and LHCb Col- 
laborations [41] at the LHC, finding good overall agreement with the CJ PDFs, but weak 
sensitivity to large- a; PDFs. 
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The large-x PDF uncertainties also affect the production rates of heavy W and Z' bosons, 
and the impact of these was studied as a function of the boson mass for Standard Model 
couplings. At high rapidity, the W'^ cross section in pp and pp collisions is mostly sensitive 
to the u quark PDF at large x and therefore well constrained. The Z' cross section displays 
a mild sensitivity to the d quark PDF at large x, reaching upwards of 40% uncertainty at 
the kinematic rapidity limit. The effects are even more pronounced for W'~ cross sections, 
where the uncertainties exceed 100% at large yw 

At central rapidity the Z' cross section is well constrained, with weak sensitivity to u 
and d quarks at the LHC for Mz' < 3 TeV. The W central rapidity cross sections, on the 
other hand, display sensitivity to large-x PDFs of the order 10 — 20% for Mw' ~ 2 TcV 
at the LHC (due to d and u) and 0.75 TeV at the Tevatron (due again to d quarks). The 
uncertainties at central rapidity directly propagate to the integrated cross sections, which 
show ^ 3% effects for Z' production, while somewhat larger for W production, amounting 
to < 20% for W'+ and < 10% for W' at the LHC for Mw < 3 TeV, and < 30% for W'+ 
and W'- at the Tevatron for Mw < 1 TeV. 

These considerations place important limits on the ability to accurately measure heavy W 
and Z' cross sections in hadronic collisions, particularly at large rapidities and boson masses 
near the kinematic thresholds of current colliders. Although our analysis is, for illustration, 
restricted to heavy vector bosons with Standard Model couplings, and the quantitative effects 
of the PDF uncertainties would be different in other models, our main point is that caution 
must be exercised when using PDFs in regions where these are not directly constrained, 
or their uncertainties underestimated, as is the case at large x. The uncertainties in the 
production cross sections can be reduced by obtaining better constraints on PDFs at large 
X, especially for the d quark. Several experiments aimed at determining the d quark PDF 
up to X 0.8 are planned at Jefferson Lab following its 12 GeV energy upgrade [21-23]. 
Uncertainties in cross sections at central rapidity, and hence in the integrated cross sections, 
will also be reduced with improved determinations of antiquark distributions at large x, 
such as the E-906/SeaQuest experiment at Fermilab [52] which plans to measure d/u up to 
X « 0.45. 

The flavor dependence of weak boson production could also be studied with pn collisions 
at the LHC or at the Relativistic Heavy- Ion Collider, with the neutron provided by a beam of 
deuterons [53]. Unlike for fixed target experiments, in a collider one can study pd collisions at 
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large positive and negative rapidities. Therefore, partons in the beam at large xi can scatter 
from partons in the target at small X2, and vice versa. In particular, measurements at large 
negative rapidity would be sentitive to quarks in the deuteron at large momentum fractions, 
offering a probe of nuclear corrections complementary to deuterium DIS, and constraining 
the nuclear uncertainties studied in this paper. 

Finally, while we have focussed on the ^/s — 7 TeV energy at which the LHC currently 
operates, in future this is planned to increase to = 14 TeV. The behavior of the cross 
sections illustrated here will not change qualitatively at the larger energy. However, for phys- 
ical W and Z bosons, the region in rapidity where sensitivity to nuclear models is greatest 
will be shifted outside of the acceptance of current experiments, limiting the usefulness of 
14 TeV data for large- a; PDF studies. On the other hand, the higher center of mass energy 
will increase the accessible W, Z' mass range (up to M^i^z' ~ 6 TeV) over a larger range 
of rapidities. 
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